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Abstract. The critical exponent <j) c , derived from the anomalous dimension of the bilinear 
operator responsible for crossover behaviour in O(N) </> 4 theory, is calculated at 0(1/N 2 ) in a 
large N expansion in arbitrary space-time dimension d = 4 — 2e. Its e expansion agrees with 
the known 0(e 4 ) perturbative expansion and new information on the structure of the five loop 
exponent is provided. Estimates of <f> c and the related crossover exponents C and 7 C , using 
Pade-Borel resummation, are provided for a range of N in three dimensions. 
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Renormalization group techniques have been widely used to study the critical properties of 
the scalar quantum field theories underlying a variety of condensed matter systems. For instance, 
Wilson and Fisher, ^, g], introduced the technique of extracting numerical estimates of critical 
exponents from the evaluation of the perturbatively calculated renormalization group functions 
at several loop orders in ^ 4 theories. Subsequently various authors have developed this method to 
very high loop orders either in fixed, (three), spacetime dimensions or in d = 4 — 2e dimensions. 
The results for various renormalization group functions at, respectively, six and five loops, which 
represents the highest orders computed, are given in f||]. In the latter case the results have 
been extrapolated to three dimensions using resummation techniques, [|, ^| fjj. These exponents 
derived by the renormalization group method are competitive with other approaches such as 
the high temperature series expansion and Monte Carlo results and are in good agreement with 
experiment. A recent and comprehensive review of the application of the renormalization group 
in this area is given in ||. Recently the critical exponent corresponding to crossover behaviour 
in 0{N) (j) 4 theory has been calculated to a new degree of accuracy in fixed dimension, |J, 
where the relevant Feynman diagrams were calculated to six loops. One motivation for that 
study rests in the realisation that the N = 5 theory of superconductivity has been observed 
in nature, pcR . Prior to this the same exponent had been computed to four loops in MS in 
d = 4 — 2e perturbation theory in [11|, which built on the lower loop calculations of [||, [Ej] . 
The resulting numerical estimate for that crossover exponent was in agreement with the high 
temperature series of [13|. One other field theoretic technique which is used in estimating 
critical exponents is the large N method where the exponents are computed order by order in 
powers of 1/N. Indeed exploiting the conformal properties of the d-dimensional Wilson-Fisher 
fixed point the technique has successfully produced the critical exponent 7] at 0(l/iV 3 ) in d- 
dimensions, [jl4j|, through use of the conformal bootstrap programme. Moreover, this method 
had developed out of the earlier d-dimensional critical point technique of [jOJ, 16] which was 
based on analysing Schwinger Dyson equations at criticality in large N. In essence the method 
efficiently reproduces the bubble summation which is the main property of 1/N expansions but 
more importantly goes beyond the leading order which the conventional bubble sum calculations 
fail to handle easily. Moreover, since the exponents are expressed as a function of d they can 
be expanded in powers of e and the coefficients compared with those of the same exponents 
computed in conventional perturbation theory. Due to the critical renormalization group the 
coefficients must be in agreement. Therefore, the information contained in the large N exponents 
can be exploited, for instance, to gain insight into the large order structure of the renormalization 
group functions at several orders in 1/N. Given the recent interest in the crossover exponent we 
will focus in this letter on its evaluation at a new order in 1/N in O(N) 4 theory. Previously 
the exponent had been calculated at 0(1/N) in d-dimensions in |L7], [l8|. To achieve this we 
follow the extension of the large N fixed point Schwinger Dyson approach of [15, 16], to the 
computation of the anomalous dimensions of composite operators, [19|. 

We recall the essential points of our calculation. The crossover exponent we are mainly 
interested in, <f) c , is computed from the anomalous dimension of the bilinear traceless symmetric 
tensor 

zab 

T ab = <rf - — <j> c <t> c (1) 
where (f> a is the field of the O(N) 4 theory and 1 < a < N, through the scaling law 

<t> c = (2 - T] c ) v . (2) 



This composite operator is relevant for a variety of critical phenomena, M, [L3L 20, Ell]. The 



exponent v has been computed in d-dimensions at 0(1/N ) in [T(|. To clarify with other work, 
H, the exponent r\ c is related to two other exponents r\ and rjo by 

r] c = V + 7] (3) 
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where r\ is the anomalous dimension of the field 4> a and has been computed at 0(1/N 3 ) in 
The remaining exponent tjq is the anomalous dimension of the bare composite operator T ab 
itself. We have chosen to express the relation for <j) c in this way since in a gauge theory the 
combination rj c would be independent of a covariant gauge parameter although the analogous r] 
and r\o would each depend on the choice of gauge. In the large N critical point method of 
the exponent tjq is extracted by inserting the operator T ab in a two point Green's function and 
extracting the residue of the simple pole with respect to the large N regularization in a well 



defined fashion according to [15, The residues of the simple pole of each Feynman diagram 
are then combined to obtain t]q. Before recalling how this regularization is introduced we note 
that the Lagrangian used in the large ./V technique is 

L = + ~ ^ (4) 

where g is the coupling constant and the field a is auxiliary. Its elimination produces the usual 
• interaction. The method of [15, [l6|, elegantly exploits the properties of the d-dimensional 



Wilson-Fisher fixed point in that, for example, the (massless) propagators of the fields of (|j) 
have simple power law behaviour. In momentum space, representing the propagator by the same 
letter as the field, the leading asymptotic scaling forms in the critical region are 

A B 

^) ~ ( k 2y- a > < k ) ~ {k 2y-p ( 5 ) 

where d = 2/i and A and B are the momentum independent amplitudes which always appear 
in the combination z = A 2 B in the computation of the Feynman diagrams. The powers of the 
propagators are related to the usual critical exponents by 

a = /j,-l + -r),(3 = 2- r}-x (6) 

where x is the anomalous dimension of the <7</> 2 vertex. It can also be determined from a scaling 
law involving v 

X = - ~ V ~ 2(/i - 1) (7) 
v 

where v is proportional to the critical slope of the coupling constant of the 0(N) nonlinear 
a model which is in the same universality class as (|j) in 2 < d < 4. In the large N critical 
point technique the propagators in the Feynman diagrams of a Green's function are represented 
by (|5|). However, in their present form when they are used to compute rjo the leading order 
large N graphs diverge. To regularize these infinities the regulator A is introduced by setting 
X — > X + A. Consequently the Feynman diagrams involve poles in A analogous to those in 
conventional perturbation theory where e in d = 4 — 2e is the (dimensional) regularization. It is 
the residues of these simple poles in A which are then used to extract tjq. It is worth stressing 
that we will compute the exponent in d-dimensions where d is arbitrary and e is not used as a 
regularization. 

For the large N renormalization of the composite operator T ab it turns out that only those 
Feynman diagrams where the operator is not within a closed <p a field loop will contribute to 
r\Q. This is a consequence of the traceless nature of the operator. Diagrams where T ab is inside 
a closed (j) a loop vanish when one computes the group theory factor of the graph. Therefore, 
at leading order, 0(1/N), there is only one Feynman diagram to calculate and applying the 



method of [19], we find 



^ = " (^2) (8) 
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where 



with 



and 
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4T(2/x - 2) 



r( M + i)r( M -i)r( M -2)r(2- M ) 

oo 

»7 = 5^ 
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Consequently, we have 



+ 



2/x77i 



+ O 



1 



(9) 
(10) 

(11) 
(12) 



which is in exact agreement with pT) , [Tq1 , though extracted with a minimal amount of effort. In 
three dimensions, (O) gives 



32 /J_ 

7T 2 iV + liV 2 



or using a Pade approximant 



1 + 



16 



(13) 



(14) 



Interestingly evaluating ( |14| ) for N = 2, 3, 5, and 16 we find the respective values 4> c = 1.105, 
1.298, 1.510, and 1.816. These are relatively close to the values obtained by other methods, |J], 
which are given in Table 1. Indeed the estimate for N = 3 is remarkably good. By contrast 
the direct evaluation of ( |i~3| ) gives respectively 0.379, 0.919, 1.352 and 1.797 indicating its poor 
convergence for low N. 



N 


4>c 




7c 


2 


1.184(12) 


0.830(12) 


0.354(25) 


3 


1.271(21) 


0.863(21) 


0.41(4) 


4 


1.35(4) 


0.90(4) 


0.45(8) 


5 


1.40(4) 


0.90(4) 


0.50(8) 


8 


1.55(4) 


0.94(4) 


0.61(8) 


16 


1.75(6) 


0.98(6) 


0.77(12) 



Table 1. Values of crossover critical exponents from 



To determine rjo,2 we have repeated the method on the 0(1/N 2 ) diagrams. Due to the way 
the large N expansion orders this would ordinarily mean that graphs up five loops would have 
to be calculated. However, when the group theory factor is computed only six diagrams remain 
with a non-zero coefficient. These are comprised of four 2-loop and two 3-loop graphs. As a 
check on our method of calculation we have redetermined v 2 from the evaluation of the exponent 
Xi using the same computer programme written in the symbolic manipulation language Form, 
[22]. The method of extracting \2 is the same as that for r]Q^ since the Feynman diagrams for 
the latter are equivalent to those for the former when the operator insertion is replaced by the 
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acf) 2 vertex. Therefore, the result of our calculation is 



VO,2 



2ix + 5 + 



14 



+ 2/i + 2 + 



(/i-2) 
1 



+ 



0* 



(/i-2) (/i-2) 2 ( M -2)3 2(/i-l) 



(15) 



where 



i?! = V'(m-i) - V>'U) 

i? 2 = V'(2/" " 3) - V'(2 - M) - ^'G" - 1) + ^'(1) 

R 3 = ^(2/i-3) + V(2-//) - V>(M-1) - ^(1) 

u' = V(2/u-2) + ^>(2-/i) - fp(fi-2) - V(2) 



(16) 



and %[)(x) is defined by VK^) = T'(x) /T(x) where T(x) is the Euler gamma function. Consequently, 
1 2//771 



(/"-I) + (/i-l)(M-2)iV 
~3,u 2 (8/i 2 -21/i + 14)i?i 



+ 



2(/i-l)(/i-2)3 



^-3) 2 
(M-l)(M-2) 3 



i? 3 + i?2 



//(V - 14/. 2 + 10m + 1) , 9 

+ ; 77777 — v - 2 + 
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41 



(M-l) 2 0x-2) 2 
4 15 



+ 



(/i-2) 
3 



(/i-2) 3 (/i-1) ( m -1)2 20* -1)3 



(/i-2) 2 
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AT2 



^ + 0( w ) (17) 



where, for completeness, we note that the values of the other exponents used to determine 

2/i(/i-l)\ , , (// 2 + /i-l) 2/i^-l) , m(3-A*)~ 



are, [15], 
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and, 1 16], 
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(/i-2) 

(2/i - 1)77! 



2/i(/i - 1) 



(/i-2) 2(/x-2) 2 



(18) 



2(/x-l) 2(/i- l)(p-2)N 

2u(2u-3) 2 



6i?i 



(/i-2) 

2(4// - 12/i 3 + 5/i 2 + 6/i - 2) , 



i?2 



+ 



m(m-i) 

3 



t/ + 4// z - 2/i + 34 + 



(/i-2) 



+ 



(/i-1) (/i-1) 2 /i 2 (/i-l) 



4(/i- l)(/x-2) 2 iV 2 



+ O 



1 

AT3 



(19) 



To check the correctness of (|l7]) we have evaluated <f> c at 0(e 4 ) in d = 4 — 2e and compared 
with the previous dimensionally regularized four loop MS perturbative calculation of the same 
critical exponent. The result (|l7]) is in exact agreement which is a non-trivial check on our 
computation since only three loop graphs are present at 0(1/N 2 ). With (17) we can expand to 
a new order in e and find 

6 C = 1 + e + e 2 + e 3 + e 4 + e 5 



8e - 8e 3 - 16(1 - C(3))e 4 - 24(1 - C(4))e £ 



N 



5 



+ 



64e - 124e 2 - 4(43 + 60C(3))e 3 + (640C(5) - 360C(4) + 976C(3) - 155)e 4 

, 1 / e 6 

+ 2(800C(6)- 1840C(5) + 732C(4) + 128C 2 (3)-144C(3) + 61)e 5 — + 0| — 



iv 2 ^iv 3 y 

(20) 

where C( n ) is the Riemann zeta function and the order symbol represents independently higher 
order terms in e and 1/N. The 0(e 5 ) coefficients will be important in future explicit five loop 
MS perturbative calculations. 

We are now in a position to examine the critical exponents in three dimensions. For the 
various ones we are interested in we have 

_ _j32_ _ 64[9tt 2 + 16] , f\_ 
9c tt 2 N 9vr 4 iV 2 \N* 

32 512 / 1 

^ ~ 3^iV " 27vr 4 iV 2 + viV 3 

•>° = ^ + (^)- < 21 > 

For reference, the other intermediate exponents are 

__8 512 / J_ 

77 ~ 3vr 2 iV 27vr 4 iV 2 + ViV 3 
32 32[27tt 2 + 1041 

" " 1 - ~ M +0 {W) ■ (22> 

In addition we record that the values of two related crossover critical exponents are 

8=1- ^ + 8] + O ( — 

32 32[9vr 2 -401 1 \ 

which are defined through the hyper scaling laws 

B c = 2fiu - <f> c , lc = 2ct> c - 2^u . (24) 

Clearly the 0(1/N 2 ) correction to (fi c is large and the series appears to diverge. By contrast 
the 0(1/N 2 ) correction to tjq vanishes in three dimensions. We have repeated our earlier Pade 
approach for <p c to see if the convergence is improved but this does not lead to a small change to 
the previous values for the exponents. This is in part due to the fact that the exponents r\ and v 
do not lend themselves to improvement by this approach. Instead one way of gaining estimates 
from our large iV results is to use the accepted values of rj and v and our value for r\o . Indeed 
in EJ the four loop estimate for <f> c was determined in an analogous fashion. Therefore, taking 
rj to be 0.033 and 0.033 and v to be 0.669 and 0.705 for N = 2 and 3 respectively, 0], we find 
the values for cp c are 1.044 and 1.196. These are in poor agreement with the respective results 
of ||. For the exponents (3 C and 7 C the large N corrections are also large for small N and each 
series appears to converge slowly. To appreciate this we have evaluated the above expressions 
for larger values of N. For N = 8 we find (j) c = 1.475, 8 C = 0.908 and 7 C = 0.567. By contrast 
when TV = 16 our expressions give 1.767, 0.977 and 0.790 for the same respective exponents 
which, by contrast, compare much more favourably with the respective values of 1.75(6), 0.98(6) 
and 0.77(12) of g]. 
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In order to improve the convergence of the series we have also examined the Pade-Borel 
resummation of the large ./V series. This involves determining the Borel function of the series 
which is defined by 



oo i ™ oo , r 

J2 a n x n = - dt e -*/« £ 
£o xJo n] 



(25) 



and then taking a Pade approximant of the integrand given that only several terms in the series 
are known. Therefore, for <j) c its Pade-Borel estimate is 



2N 



-Nt 



dt 



[1 - ait + (a\ - \a 2 )t 2 } 



where 



Ol 



7T 



32[9vr 2 + 16] 
9^ 



(26) 



(27) 



We have evaluated the integral numerically for various values of N and recorded the results 
in Table 2 where the estimates for (3 C and 7 C by the same method are also given. The final 
column is the sum of the estimates in the second and third columns and represents another way 
of estimating C through the scaling relation since we have noted that the large N series for </> c 
appears to diverge rapidly for low N. For N > 4 the large N estimates for <j) c and the sum 
Pc + 7c are in fairly reasonable agreement. For N = 2 and 3 the estimates undershoot those 
of [H| though the combination (3 C + 7 C is closer. For the other exponents the values for (3 C are 
competitive for N > 5 whilst those for j c appear to be in good agreement for the lower range 
of N. 



N 


4>c 




7c 


/3 C + 7c 


2 


0.988 


0.664 


0.367 


1.031 


3 


1.187 


0.768 


0.459 


1.227 


4 


1.323 


0.830 


0.529 


1.359 


5 


1.422 


0.871 


0.582 


1.453 


8 


1.603 


0.934 


0.689 


1.623 


16 


1.790 


0.980 


0.817 


1.797 



Table 2. Pade-Borel estimates of crossover exponents. 



In conclusion, we have provided the 0(1/N 2 ) corrections to a set of crossover exponents 
related to the composite operator T ab in O(N) </> 4 theory. Although the leading order exponents 
could be summed to give numerical estimates which are competitive with explicit perturbative 
calculations in three dimensions the new higher order correction indicate that the series are slowly 
converging. Applying the Pade-Borel resummation technique generally improves the estimates 
in comparison with the results of || though it ought to be borne in mind that 0(1/N 2 ) results 
represent only three terms of a series in contrast to || which analysed six terms of a series. 
Nevertheless since the critical exponents are computed in d-dimensions they will complement 
future higher order perturbative calculations and, further, the large N method can equally be 
applied to the determination of crossover exponents of bilinear and other composite operators 
to the same large order in this and other scalar quantum field theories which underpin critical 
phenomena. 
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